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Abstract 

We investigate homogeneous and isotropic cosmological solutions sup- 
ported by the SU(2) gauge field governed by the Born-Infeld lagrangian. In 
the framework of the Friedmann-Robertson- Walker cosmology with or with- 
out the cosmological constant A, we derive dynamical systems that give a 
rather complete description of the space of solutions. For A = the effec- 
tive equation of state s{p) is shown to interpolate between p = — e/3 in the 
regime of the strong field and p = e/3 for the weak field. Correspondingly, the 
Universe starts with zero acceleration and gradually enters the decelerating 
regime, asymptotically approaching the Tolman solution. 
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1 Introduction 



It is tacitly assumed that the large scale massless Yang-Mills (YM) fields, which 
could exist in the early Universe before phase transitions, do not play any signifi- 
cant role in cosmology. Partially this is related to the scale invariance of the YM 
lagrangian implying that primary YM excitations would be diluted during infia- 
tion. String theory suggests the Born-Infeld (Bl) type modification of the YM 
action which breaks the scale invariance thus canceling the above objection. There- 
fore, it seems reasonable to investigate the YM cosmology with the non-Abelian 
Born-Infeld action. In the existing literature one finds several papers discussing 
cosmological models with the U(l) BI matter 0, H. Such models are necessarily 
anisotropic (or inhomogeneous) , since there is no homogeneous and isotropic con- 
figurations of the classical U(l) field. Non-Abelian BI (NBI) cosmology was not 
studied so far. Here we investigate this problem in the framework of the traditional 
Friedmann-Robertson- Walker (FRW) cosmology. 

A notable property of the SU(2) Yang-Mills field is that it admits homogeneous 
and isotropic (Ge-iiivariant) configurations. Indeed, the energy can be distributed 
between three color components in such a way that the resulting stress-tensor is 
compatible with the maximal symmetry of the three-dimensional space. In the case 
of S^, the corresponding ansatz was given in ||^, |^, |^, two other cases (hyperbolic and 
flat) were treated in . An interesting feature of these configurations is that they 
produce a varying with time Chern-Simons density, which might lead to topological 
fermion non- conservation P, 0. As was shown in [0, in the case of the spatial 
geometry, the YM field generically oscillates between two neighboring topological 
sectors of the YM theory, while the limiting unstable static solution plays a role 
of the cosmological sphaleron p|, Here we do not discuss this issue further 

but rather concentrate on the aspects related to the scale invariance breaking in 
the NBI theory. An immediate consequence is that the equation of state p = e/3 
arising in the usual YM theory |^ (mimicking the photon gas) changes to some 
more complicated equation which now allows for a negative pressure. One could 
wonder whether an accelerating expansion becomes possible, that is, whether the 
sum e + 3p can be negative. It turns out that the lower limit of acceleration which 
can be achieved within the present model is precisely zero: when YM field strength 
is much larger then the BI 'critical field', the equation of state is p = —e/3. Such a 
state equation is typical for an averaged distribution of strings, where its origin is 
quite simple: each string has a one- dimensional equation of state p = —e, averaging 
over all directions in the three-space gives p = —e/3. The deeper reason for this 
coincidence is not fully clear for the moment, though apparently it is related to the 
fact that the origin of the BI lagrangian lies in the string theory. 

As it was widely discussed recently JIT], |T2|, |T3|, |T?|, the definition of the NBI 
action is ambiguous. One can start with the U(l) BI action presented either in the 
determinant form ^ 

S = TT— J \l- det(5f,,„ - F„„) d'^x 



167r 
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or in the 'resolved square root' form 




167r7 V 2 16 ^ ^ 

In the non-Abehan case, the trace over color indices must be specified. One partic- 



ular definition, suggested by Tseytlin |12|, is a symmetrized trace. It prescribes a 



symmetrisation over all permutations of the gauge matrices in an expansion of the 
determinant in powers of the field strength before the trace is computed. Inside the 
symmetrized series expansion, the gauge generators effectively commute, so both the 
determinant and the square root forms are equivalent (which property may serve 
as an alternative definition of the trace operation). This property does not hold for 
other trace prescriptions, e.g., an ordinary trace. In this latter case, it is common 
to apply the trace to the square root form. Note, that the string theory favors the 
symmetrized trace definition once lower orders of the perturbation theory are taken 
into account |T2|, [T^, while higher order corrections seem to violate this prescrip- 



tion |T6|, [T^, 0, |T9|- Here we choose the 'square root /ordinary trace' lagrangian just 
for its simplicity. It is worth noting that in other (somehow related) problem of 
sphalerons in the NBI theory, discussed recently both in the ordinary pOl and sym- 



metrized trace [^Tl versions, basic qualitative features of the solutions turned out 
to be the same. We realize, however, that the distinction between different trace 
prescriptions may become significant near the cosmological singularity. 

Fortunately, the dynamics of the scale factor in the SU(2) Einstein-NBI theory 
(with ordinary trace) decouples from the whole system of equations and can be 
fully analyzed using the dynamical systems approach. We perform this analysis for 
closed, spatially fiat and open FRW models both for zero and non-zero cosmological 
constants. Given the solution for the scale factor, the field equations for the YM 
field can be integrated in terms of elliptic functions. 



2 SU{2) Ge-invariant connection 

We start with the action 

5 = -^/{^^ + /5'(7^-l)} V^d\ (1) 
where R is the scalar curvature, (3 is the BI critical field strength, and the quantity 

corresponds to the 'square root/ordinary trace' NBI lagrangian. Throughout the 
paper we consider cosmological models with the homogeneous and isotropic three- 
space, which will be presented as usual: 
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where E is one of the functions sinr, r or sinhr depending on the type of the 
spacetime: closed, spatially flat or open (labeled as = 1,0, —1 correspondingly). 

In what follows, we deal with the homogeneous and isotropic configurations of 
the SU(2) YM field. Although the construction of the SU(2) YM ansatz for all 
FRW models was discussed earlier |^, it is worth to reconsider it, since in [0] the 
conformal invariance of the YM lagrangian was explicitly used. In fact, the result 
obtained in |^ was just the Gg-invariant connection independent of the choice of 
the action for the gauge field. Here we rederive the Gq ansatz without any appeal 
to the scale invariance. 

The homogeneous and isotropic ansatz for the YM field can be obtained starting 



with the Witten spherically symmetric ansatz p2| and then extending the symmetry 
to the full Gq group. The Witten ansatz reads: 

A^dx^" = UoTr dt + uiTr dr + [Ka^e - (1 - Ki)T^] d9 

+ [(1 - K^)Te + K2T^] sin 6 d<P, (3) 

where four functions ojq^uii, Ki, and K2 depend on time t and the radial coordinate r 
(spherical coordinates r,9,4) are understood), and the coordinate dependent basis 
in the color space is introduced as follows 

Tr = Ti sin 9 cos + T2 sin 6* sin + T3 cos 9, Tg = dgTr, sin 9 = d^jjTr. (4) 

These generators satisfy the standard normalisation condition and the commutation 
relations 

tr TaTt = hab, [Ta, T,] = te'^^'^n. (5) 

From four functions entering the ansatz only three are physical, while a{r,t) = 
aicta.n{K2/ Ki) can be gauged away by a gauge transformation preserving the 5*0(3) 
symmetry. In the following it is convenient to choose a parametrisation 

Ki = / cos a, K2 = / sin a, 

= Qq + a, Ui = Qi + a'. (6) 

where a prime and a dot denote derivatives with respect to the radial variable and 
time. 

In order to extend the symmetry to Gq, we consider the following gauge invariant 
tensor: 

= f;,f/:\ (7) 

Since in the homogeneous and isotropic three-space the only mixed second rank 
tensor is the Kronecker delta, one should have for the spatial part of P^''^: 



Pi = Pit) 6: 



Using (0-|^, we obtain the system of equations for f2o,i, /, and a, which is solved for 
the first three functions in terms of a single function of time w{t), while a remains 
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arbitrary. The following gauge leads to a maximal simplification of the field strength: 

= - "^y.^."^^! , (8) 



1 + S2(w;2 


-k) 




-k) 


1 + S2(w;2 


-k) 



a = arctan 
In this gauge, the field tensor reads 



/=^1 + S2(^2_^)^ 



VI - kT? 



(9) 
(10) 

(11) 



T = w{Trdt A dr + TeJ:dt A d9 + T^E sine dt A d^) 

+J:{w^ - k) (T^ drAdO- Tg sin 9dr Ad^ + sin 6 dO A d^) . (12) 

3 FRW cosmology 

Consider the FRW cosmological models parameterizing the spacetime metric as 

ds^ = NHe -a'df. 

Substituting (p!2D into Eq. (|I]), we obtain the following one- dimensional action: 

" 3 a(fcAr2 _ ^2) 



S^ = — I dt 
Arc 



2G 



N 



Na'(3\n 



where now 



7^ 



1 



3^2 



3{w^ - ky 9u;2(w2 - ky 



The lagrangian contains two dimensional parameters: the Newton constant G 
and the BI 'critical field' p. After a coordinate rescaling t P'^^'^t, a — > P~^^'^a, 
which makes all quantities dimensionless, the reduced action becomes 



^1 



1 f , (3a(kN^ 



N 



-gNa^ 



where 



K 



and V 



(1 -K2)(l + y2) 
k) 



(13) 



aN 

and g = pC is the remaining dimensionless coupling constant. 

Variation of the action over gives a constraint equation; after obtaining it, we 
fix the gauge = 1. In this gauge the constraint equation reads 



ga' {V-l)--{d' + k) 



0, 



(14) 
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where 

This expression can be rewritten in the standard form 

k SirG ,^ , 

e, (16) 



c? 3 

where the energy density is given by 

e = e,{V-\), (17) 

with ec = (3/4:71 playing a role of the BI critical energy density. 

Variation of the action with respect to a gives the acceleration equation 

Again, rewriting it in the standard form (using (|1^)) 

d AirG , . , , 

- = ^(^ + 3p), (19) 

a 6 



we can read off the pressure 



{S-V- 2V-') . (20) 



^ 3 

Now comparing (|1^ and (|2y) we obtain the following equation of state 

3 [Sc + e) 

It is worth noting that the BI critical energy density corresponds to the vanishing 
pressure. For larger energies the pressure becomes negative, its limiting value is 
p = —e/S. In the opposite limit e <^ Ec one recovers the hot matter equation of 
state p = e/S, reflecting the scale invariance of the YM action to which the NBI 
action reduces at low energies. 

Finally, the variation over w gives the YM (NBI) equation 

.. w(k - w^) f 1 - 2aw (I 1-K^\ 



a 



From this one can derive the following evolution equation for the energy density: 

e = -2^ £i£±M, (23) 
a e + Ec 

which can be easily integrated to give 

a^{e + 2ec)e = const. (24) 

From this relation one can see that the behavior of the NBI field interpolates between 
two patterns: 1) for large energy densities {e S> Ec) the energy density scales as 
6 ~ a"^: 2) for small densities e <^ Ec one has a radiation law e ~ a~^. 
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4 Spacetime evolution 



Using the constraint equation (|T^), one can express e in terms of tlie scale factor a 
and its derivative. Substituting it into Eq. (p!8|) we obtain a decoupled equation 
governing the spacetime evolution: 



2ga{a^ + k) 
' 2ga? + 2,{a? + k)' 



(25) 

The condition of positivity of the energy density in (JT^) defines a boundary in the 
phase space: 

> -k. (26) 

Let us analyze the evolution equation (^) by means of the dynamical systems 
tools. For this, it is convenient to introduce h = a and pass to another independent 
variable r, such that dt = hdr, where h = 2go? + 3(6^ + k). This leads to the 
following dynamical system: 

a' = hh, b' = -2ga{b'^ + k), (27) 



where primes denote derivatives with respect to r. Notice that the system ( ^71) is 
invariant under a reflection (a — —a, b —b), so further, without loss of generality, 
we discuss the expanding solutions. 

First of all we establish that the dynamical system ( pTj) admits a first integral. 
Considering a as a function of b, one can write: 

da hb 



db 2ag{lP + k)' 
This equation can be easily integrated to give 

3 (6^ + A;) % Aga^ (b^ + k) = C, (28) 

with some constant C. From this relation one can see that all solutions inside the 
physical boundary (^6]) cross the hne a = 0, so that the singularity is unavoidable. 
Another observation is that, unlike the standard hot FRW cosmology, d remains 
finite while a tends to zero. 



Though the integral (|28|) allows one to draw the complete phase portraits of (|27D , 
let us take a look at the stationary points of the dynamical system. These are 
different for different k: 

• k = 1, closed universe. The only singular point is a = 0, 6 = which is a center 
with the eigenvalues ±i^/6g, see the phase portrait in Fig. All solutions 
are of an oscillating type: they start at the singularity (a = 0) and after a 
stage of expansion shrink to another singularity. It is easy to see from 
that the functions a and d remain bounded for all solutions. 



^In this and the following figures, solutions evolve from left to right in the upper half-plane as 
time changes from — oo to oo, and from right to left in the lower half-plane. 
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• = 0, spatially flat universe. There is a singular line 6 = each point of which 
represents a solution for an empty space (Minkowski spacetime). This set is 
degenerate, and there are no solutions that reach this curve for finite values 
of a (see Fig. || for the phase portrait). All solutions in the upper half-plane 
after initial singularity expand infinitely. 

A remarkable fact is that for this case one can write an exact solution of (^Tf ) 
(in an implicit form): 

4:y/g{t -to) = V3(n- arctanr]"^ + vr/2) , 

where = ^/2a/ \J \/a^ + C — a"^. The metric singularity is reached at t = to- 

• k = —1, open universe. There is a center at a = 0, 6 = with the eigenvalues 
±iy/6g, but it lies outside the boundary of the physical region. Other singular 
points are (a = 0, b = ±1) (Fig. H). These points are degenerate and cannot 
be reached from any point lying in the physically allowed domain of the phase 
plane. The only solutions which start from them are the separatrices b = ±1 
that represent (part of) the fiat Minkowski spacetime in special coordinates. 

One can easily see that all solutions in the upper part of the physical domain 
a > 1 start from the singularity and then move to a ^ oo, a — 1. 

The global qualitative behavior of solutions does not differ substantially from 
that in the conformally invariant YM field model, except near the singularity. One 
can find the following power series expansion near the singularity a(0) = 0: 

where 6o is a free parameter. Absence of the quadratic term means that the Universe 
starts with zero acceleration. This is what can be expected in view of the equation 
of state p ~ —e/3 at high densities. 

For large a, the dynamics of the system (|18]) approaches that of the hot FRW 
models (for small energy densities one recovers the equation of state of radiation). 

5 Cosmological constant 

Here we extend our analysis to the case of non-zero cosmological constant A. This 
leads to the following one- dimensional action: 

^ {kN^ - a^) + gNa' (l - ^(1 - K^){1 + V-) - . 

Clearly, the equations governing the gauge field dynamics (^) are unaffected by 
the cosmological term. The metric equations can be decoupled again using the 



^1 



AT^Gf5 



dt 
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constraint equation (jM]): 

The physical region of the phase space is defined by positivity of the energy of the 
gauge field (0), now the boundary being 

-\a^ + k> 0. (31) 

For the closed and spatially fiat universes and a negative cosmological constant 
the allowed domain coincides with the whole plane a, a. In other cases, there is a 
boundary corresponding to the solution with zero gauge field, i.e., the (anti) de-Sitter 
space. 

The cosmological constant gives rise to new types of solutions including de-Sitter- 
like, which are non-singular. However, in the case when there is a singularity, its 
structure is determined by the leading term in the equation of state, namely, by 
the BI pressure, and is thus unaffected by the cosmological constant. The generic 
solution near the singularity satisfies the following series expansion 

a(t) =M-^^i^feot' + 0(t^), t^O, (32) 

io 

where 6o is a free parameter. Again, the Universe starts with zero acceleration. 

To proceed further with the analysis of Eq. (|5D[), we rewrite it as the following 
two-dimensional dynamical system: 



a' 



{2g - X)a^ + 3(6^ + k)] b, b' = \X{Ag - X)a^ - 3(6^ + k){2g - A) 



a, 
(33) 

where primes denote derivatives with respect to a new time variable r defined by 
dt = 3[{2g — X)a^ + 3(6^ -|- fc)](ir. The system is symmetric under a simultaneous 
refiection of both variables, so we can think about a half of the phase space. 

Similarly to the previous section, one may treat 6 as a function of a and thus 
obtain the first integral, which completely defines the integral curves of (|33|) : 

+ Aga^ {b'' - ^a" + k^=C. (34) 

Note that the terms in brackets are non-negative within the allowed domain (|3T|) . 

The structure of the solution space can be revealed using this first integral. Let 
us start describing singular points of the dynamical system for different k. 

• k = 1, closed universe. For a negative cosmological constant the only singular 




point is (a = 0, 6 = 0), which is a center with eigenvalues ^ = ±2>i^?>{2g — A). 
All phase space trajectories are deformed circles, and all solutions start from 
the singularity and reach the singularity in the future. 
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For a positive cosmological constant, the allowed domain is bounded by a 
hyperbola 6^ — > —1 which represents the de-Sitter space. 

For relatively small A the point (a = 0, 6 = 0) remains a center. At the same 
time, another pair of singular points appears: 



3(2^ - ^) ^ _ (35) 



For \ < 2g they lie within the allowed domain and correspond to the static 
Einstein universe. These are saddle points with eigenvalues 



QVQg I X{2g-X) 
^ = ^^V Ag-X ■ 

Entering them separatrices correspond to either a solution developing from 
the singularity into the static universe, or one rolling down from an infinite 
radius to the static universe (Fig. |). These separatrices divide the phase 
space into the domains containing different types of generic solutions. Near 
the origin, all solutions evolve from an initial to a final singularity. In the 
region near the physical boundary, the solutions are non-singular and evolve 
for an infinite time first shrinking to some finite value of the scale factor and 
then ever expanding. Generic solutions of the third type possess a singularity 
but are non-periodic: the universe expands forever. 

With A further increasing, the saddle point approaches the origin and finally, 
when X > 2g, swallows it. The character of the singular point (a = 0, 6 = 0) 
changes — it becomes a saddle point ((Fig. ||). The separatrices entering it 
are solutions of an infiationary type. They divide all generic solutions into 
two classes: nonsingular de-Sitter-like (i.e., ever-expanding in the future and 
ever-contracting in the past) and ever-expanding solutions with an initial sin- 
gularity. When X > Ag, the Eq. (^) has a real singular point lying outside 
the physical region. 

k = 0, spatially fiat universe. The only singular point is the origin {a = 0,b = 
0) which now is degenerate. For A < the allowed domain is the whole plane. 
All solutions are of an oscillatory type evolving from the initial to the final 
singularity for a finite time. 

For a positive cosmological constant the allowed domain consists of the upper 
and lower parts of a cone fe^ — ^a^ > whose boundary corresponds to the 
(part of) de-Sitter space described in the infiating coordinates. All solutions 
that lie within the allowed domain have an initial singularity and are ever 
expanding in the future. 

k = —1, open universe. The allowed domain ([31| ) lies outside of either an 
ellipse (A < 0), or a hyperbola (A > 0). 
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In this case two types of singular points — the origin and 



a = ± 



3{\-2g) 
^ A(4^7-A)' 







(36) 



(when they are real) lie outside the allowed domain and hence are of no interest. 

Another pair of singular points is (a = 0, 6 = ±1). These are degenerate (as 
in case with zero cosmological constant). Since they lie on the boundary of 
the allowed domain, the only physical solutions approaching them correspond 
to zero energy of the BI field ((anti)-deSitter). 

All generic solutions within the allowed domain possess a singularity and are 
either oscillating (for A < 0) or ever expanding (for A > 0). 



6 Dynamics of the gauge field 

One can easily see that the gauge field influences the metric only through the quan- 
tity r (|15D, related to the energy density. This quantity obeys the differential 
equation: 

r = 2-(^-v], (37) 



a \V 

which follows from (^2]) and does not depend either on a detailed YM dynamics, or 
a particular gauge choice (e.g., = 1, or = a). This equation can be integrated 
once, giving V function of a: 



V^^l,s{^y, (38) 

where is an integration constant, and a numerical coefficient was introduced for 
future convenience. This equation gives us a possibility to fully describe the gauge 
field dynamics. Recalling the definition of V (plSj) , one can separate the metric and 
the gauge field variables as follows: 

(39) 



Since the right hand side of this equation is strictly positive, we can find the domain 
of variation of the function w. It is symmetric with respect to w ^ —w, except for 
the case oq < 1, = 1, when w oscillates near the value w = 1 (or w = —1) without 
crossing the axis w = 0. 

The equation (|39|) can be solved in terms of the Jacobi elliptic functions, similarly 
to the case of an ordinary YM lagrangian ES] : 



'k + 

w{t) = ^k + alcn{ V2ao{T - Tq), ^ ] , (40) 
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where tq is the integration constant, and the argument to the Jacobi function is 
defined as 

'•i a{t')N{f) dt' 



A generic solution for w{t) is of an oscillating type. In the case ao < 1, /c = 1, the 
solution, oscillating near one of the values w = ±1, does not cross the line w = 0. 
In all other cases solutions oscillate around the origin. The effective frequency of 
oscillations is determined by the rate of growth of the argument r. To compare 
the situation with the ordinary EYM cosmology, we notice that in this case the 
solution for w{t) is still given by (^OD, but with a different definition of the phase 
variable r(t) ^ 

Nit') 



a{t') 

Clearly, the main difference with the ordinary EYM cosmology relates to small values 
of a. One can see that near the singularity (a 0) the YM oscillations in the NBI 
case slow down, while in the ordinary YM cosmology the frequency remains constant 
in the conformal gauge = a, or tends to infinity in proper time gauge N = 1. 
The behavior of the YM field for closed and spatially fiat models is illustrated in 
Figs.(i^. 

When the solution starts from the metric singularity, the gauge field function 
has the following series expansion (the metric is supposed to satisfy the expansion 

m))- 

I ^0"^ f2 I ( gXMk-w^) boa{\-2g) \ , , 



where a = ±y3{k + 6q)^ — Ag'^ik — Wq)"^, Wq is a free parameter (as well as the sign 
of a), while bo is a parameter from the expansion of the metric (|32|). 



7 Discussion 



In spite of considerable complications which the Born-Infeld non-linearity intro- 
duces to the Einstein- Yang-Mills coupled equations, the dynamics of the FRW NBI 
cosmology turned out to be separable and admitting a rather complete analysis in 
terms of the dynamical systems theory. This is due to the fact that the NBI YM 
equation admits a first integral expressing the variation of the energy density under 
cosmological evolution. Thus one obtains a decoupled equation for the scale factor 
both without or with the cosmological constant. We have given a rather complete 
classification of possible solutions for closed, spatially fiat, and open universes in 
these cases. 

One of the intriguing features associated with the NBI theory is the possibility of 
a negative pressure. Namely, the pressure is negative when the field energy density 
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exceeds the 'critical' BI energy density. This is still insufficient to mimic an inflation, 
since the lowest value is only p = —e/3, which coincides with the equation of state 
of the gas of Nambu-Goto strings in three spatial dimensions. In the framework 
of the FRW cosmology, the equation of state of the NBI YM fleld continuously 
interpolates from p = —e/3 near the singularity to the 'radiation' equation p = e/3 
at large time. Correspondingly, the energy density is evolving near the singularity 
according to the law a~^. The FRW NBI universe starts with zero acceleration and 
achieves the Tolman expansion regime, when the energy density is diluted sufficiently 
to make the Born-Infeld non-linearity negligible. 

Once the scale factor evolution is determined, one is able to flnd an analytic 
solution for the YM function is terms of elliptic functions. Qualitatively the YM 
dynamics remains the same as in the usual YM theory, though we have observed 
that near the singularity the oscillations of the YM fleld are damped. 
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Figure 1: Phase portrait of the dynamical system for the closed NBI universe. 
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Figure 3: Phase portrait of the dynamical system for the open NBI universe. Phys- 
ical domain in this case is — 1 < 6 < 1. 
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Figure 4: Phase portrait of the dynamical system (^) for the closed NBI universe 
with a cosmological constant < X < 2g. Physical domains are bounded by hyper- 
bolas plotted in dashed line. 
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Figure 5: Phase portrait of the dynamical system (^) for the closed NBI universe 
with a cosmological constant A > 2g. 




Figure 6: Sample solution for the closed universe filled by the NBI YM field. Dotted 
line — gauge field function w, solid line — metric function a. The values of parameters 
are g = \, C = b, wq = 1.3. 
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Figure 7: Sample solution for the spatially flat universe with NBI YM field. Dotted 
line — gauge field function w, solid line — metric function a. 
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